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FOREWORD 

This document is an interim report describing results of the 
first ten months of effort on Contract NAS8-31440. ’‘Space Shuttle 
Contamination Due to Backflow from Control Motor Exhaust." The 
contract was initiated on 10 April 1975 as a 10 -month effort to de- 
velop analytical techniques for predicting the back-scattering of 
outgas contamination from an orbiting spacecraft. 

The Scope of Work of the contract has been modified to include 
additional tasks and the period of performance extended for an addi- 
tional nine months of effort. The final report is being deferred until 
the end of the extended period of performance. 

This study was performed by personnel of the Lockheed Missiles 
& Space Company, Inc., Huntsville Research &: Engineering Center, for 
the Space Sciences Laboratory of NASA-Marshall Space Flight Center. 
The NASA technical monitor for the study is Dr. R, J. Naumann. Dr. 
Naumann made significant contributions to this study through his own 
research efforts and through enlightening discussions with the author. 
His interest and assistance are greatly appreciated. 
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Section 1 

INTRODUCTION AND SUMMARY 

It has long been recognized that the quality and resolution of optical 
astronomical observations are limited by the disturbing effects induced 
during passage of light through the earth’s atmosphere. For this reason, 
measurements have been planned to be performed on board orbiting space- 
craft at altitudes where the disturbing influence of the earth’s atmosphere 
is negligible. Unfortunately, even in this near -vacuum environment, optical 
measurements can be degraded by the atmosphere induced by the spacecraft 
and instrumentation packages. A primary concern is the condensation and 
adsorption of outgas products on the various optical surfaces (lenses, mirrors, 
etc.) such that the optical properties of these surfaces are altered. Extremely 
small amounts of contamination can significantly affect the transmissivity and 
reflectivity of these optical surfaces, thus seriously degrading the intended 
astronomical measurements. 

Outgas products can reach sensitive optical surfaces by: (1) direct flow 
from sources of outgas contamination; (2) scattering off spacecraft surfaces; 
or (3) back- scattering due to intermolecular collisions. The first two modes 
of transport are rather easily analyzed to obtain contamination flow estimates. 
The third mode is somewhat more difficult to analyze because of the compli- 
cated collision mechanics. 

This study is addressed to the analysis of back- scattering of outgas 
contamination due to intermolecular collisions. The purpose is to develop 
analytical tools for making reasonably accurate quantitative estimates of 
outgas contamination return flux, given a knowledge of the pertinent space- 
craft and orbit conditions. Two basic collision mechanisms were considered: 

(1) collisions involving only outgas molecules (self -scattering); and (2) collisions 
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between outgas molecules and molecules in the ambient atmosphere (ambient- 
scattering). For simplicity, the geometry was idealized to a uniformly out- 
gassing sphere and to a disk oriented normal to the freestream. The method 
of solution involved an integration of an approximation of the Boltzmann kinetic 
equation known as the BGK (or Krook) model equation (Ref, 1). Results were 
obtained in the form of simple equations relating outgas return flux to space- 
craft and orbit parameters. Results are compared with previous analyses 
based on more simplistic models of the collision processes. 
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Section Z 

MODEL DEVELOPMENT 


The two basic mechanisms involved in back- scattering of outgas con- 
tamination are illustrated in Fig, 2-1. One mechanism involves collisions 
between the outward flowing outgas molecules such that some of the mole- 
cules are deflected back into the outgassing surface. The other mechanism 
involves collisions between the outgas molecules and molecules in the ambient 
atmosphere. The atmospheric molecules are moving at orbital velocities 
relative to the spacecraft and are thus highly energetic compared to the out- 
gas molecules. 


Molecular flow processes of this nature are described by the Boltzmann 
kinetic equation which describes the time and spatial rate of change in the 
velocity distribution functions as a result of intermolecular collisions. The 
full Boltzmann equation is given in Eq.(Z.l) for one of the species in a binary 
(two species) gas; 


afi 

at 


+ V 





v\ 


duO 




( 2 . 1 ) 


where the f s are velocity distribution functions for the separate species, the 
ip's are differential collision cross sections, and 





— k. — > 
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Fig, Z-1 - Return Flux MechanisTi for Outgas Contamination 
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= fl(vi), f’i=fi(v\), F^ = fj(Vj), F'j=f^lV'j) 
^2 = ^ 2 <^ 2 )' ^2 = ^ 2 <^ 2 ) 


The single subscripts "1" and "2" denote the individual species in the binary 

mixture. The double subscripts refer to collision parameters with the " 1 1" 

subscripts referring to self-collisions and the " 12" subscripts referring to 

cross collisions. The parameters in the equation include the molecular velo- 

city vector v, position vector r and time t. The prinies refer to the "other" 

molecule in self-collisions. The capital letters refer to velocities after the 

collision. The relative velocities v and v refer to relative velocities 

^11 ^12 

between the colliding molecules prior to the collision. The solid angle element 
dco contains the line of centers between the colliding molecules at the time of 
impact. 

The left-hand side of Eq.(2.1) describes the rate of change in the 
distribution function, and the right-hand side contains the collision integrals 
which account for the effects of intermolecular collisions. When the right- 
hand side is zero, the flow is free molecular, and the left-hand side can be 
solved rather easily for a variety of problems. When the right-hand side is 
not zero, however, a rigorous solution would require an evaluation of the 
collision integrals at all points in the flow field, a formidable task which 
generally requires an enormous amount of computational effort. The only 
known exact solution of the Boltzmann equation with collisions is the 
Maxwellian distribution function toward which f relaxes at equilibrium. 

In order to allow practical solutions of molecular flow processes with 
collisions, the BGK model equation (Ref. 1) was developed in which the collision 
integrals in the full Boltzmann equation were approximated by simple relaxa- 
tion terms which simulate the relaxation of the full Bcltzmann equation to the 
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Maxwellian equilibrium distribution. The BGK kinetic equations for a binary 
mixture are given as follows (Morse, Ref. 2); 


^ + (^11 - 
8 r 


9^2 - 

-gr + V • — = 1/22 ^^22 " ^2^ ‘^21 ^^21 ‘ 


where the ii subscripts indicate parameters for self-collisions, and the ij 
subscripts indicate cross-collisions. The symbol u represents collision 
rates and M represents equilibrium Maxwellian distributions defined by 



(2.3) 


By taking moments on Eqs.(2.2) for the conservation of mass, momentum 
and energy, the following relations are obtained for the model pt^rameters: 


n. 




(2.4a) 


trii 0*^2 ""2 ^21 = “l ^ ^2 ^2 

ni Uj = nj ujj =y*v fjd^ V 

n^ U2 = ^2 "22 =y u f2 V 






(2.4b) 


(2.4c) 


(2.4d) 


(2.4e) 


2-4 


LOCKHEED - HUNTSVILLE RESEARCH & ENGINEERING CENTER 



r 


I 

1 




LMSC-HREC TR D496676 


I ’'"l ®1 




m, 


f ‘‘i 


(2.41) 


2 ^"2®2 




(2.4g) 


These conservation equations alone suffice to determine the self-collision 
parameters and 0^. in addition to the densities m. Additional constraints, 
however, are required to determine the cross -collision parameters u^^ and 
0.j. Morse (Ref. 2) and Hamel (Ref. 3) required that the relaxation of momen- 
tum and temperature differences between the two species be the same as that 
derived from the full Boltzmann collision integral for Maxwell molecules. 

This leads to 


^12 " ^21 "" ^"'l ^1 ^Z "^ 2 ^ 


(2.5a) 


2m, m, 

®12 = ®1 + — [(»2-®l> ^-6^(-2--l) ] 


(2.5b) 


(mj + m^) 


2m, m, ^ 

®21 ' ®2 * *2 ['®1 -» 2 > + -bT <" 2 -“ l > 1 

(m^ + m^) 


(2.5c) 


It should be noted at this point that the developers of this kinetic model made 
use of the model for near equilibrium conditions, whereas the present appli- 
cation is for a decidedly non-equilibrium case. Willis (Ref. 4) cautioned about 
the use of the Krook model in highly non -equilibrium conditions such as these 
with his example of the calculation of sphere drag in hyperthermal, near free 
molecular flow. Before proceeding further in the use of the model and the 
model parameters defined in Eqs. (2.5), a check should be made of its validity 
in the present application. This was done by calculating directly the mean 
flow and kinetic temperature parameters for hyperthermal scattering of hard 
spheres and comparing with the relations in Eqs. (2.5). For hyperthermal 
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flow, the outgas molecules (subscript 1) may be conside: ed stationary and 
the ambient molecules (subscript 2) moving at the single velocity Equa- 

tions (2.5) then become 


m- 


^21 "" m ^ + m 2 ^*2 


u. 


12" 3k 


2 


(m^ + m2) 

ft J_ 2 

^21" 3k ■; ; T 2 


(mj +m 2 ) 


(2.6a) 


(2.6b) 


(2.6c) 


The mean flow velocities and kinetic temperatures for the hyperthermal 
scattering of hard spheres can be shown to agree exactly with the relations 
in Eqs. (2.6). These model parameters, therefore, can be used with great 
confidence in the present application. 

The remaining parameters to be determined are the collision frequencies 
The appropriate values of these parameters depend on the type of collision 
process under consideration. For the present problem, they depend on whether 
the collision process is self-collisions or ambient-collisions. This determina- 
tion will be made later as we consider each process separately. As a starting 
point, however, we can bear in mind that, for near -equilibrium flows, in which 
the distribution functions are small perturbations from the Maxwellian, taking 
moments on the single species BGK equation yields the Navier -Stokes equations 
with the correct viscosity -shear stress relation when (Ref. 5): 


V = kTn//i (2.7) 

where T and n are the local temperature and number density, respectively, 
and is the Chapman-Enskog viscosity. This value has been used by a number 
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of investigators for a wide variety oi problems, even for highly non- 
equilibrium cases. 

In our analysis of the return flow problem, we assume nearby free 
molecular flow conditions. The distribution function for the outgassing 
species is therefore a small perturbatl n from the free molecular solution: 


A 



( 2 . 8 ) 


where fj q is the fi ee molecular solution and f^ is the perturbation due to 
collisions. For velocities directed away from the oatgassing surface, fj << 
f^ Q. Since return flow is due only to collision processes, (we consider only 
convex or flat surf^.ces), the free molecular solution f ^ ^ is zero for velo'^ities 
directed toward the outgassing surface. Applying these relations to the BGK 
equation for the outgassing species (Eq.(2.2a)) yields: 


Ofj ^ 9fj 

.j V • — - 1^11 (Mjj - £j|) + 1^12 ^^J2 ■ 

dr 


(2.9) 


for mv.iecules flowing toward the se.rface. Also bat* ;d on the small perturba- 
tion assumption, the parameters involved in and are obtained by 

using the free molecular distribution functions f ^ ^ and f^ q in Eqs. (2.4). 

Note that Eq. (2.9) contains contributions to the return flow distribution 
fimction due to both self -collisions and ambient collisions, and that the sepa- 
rate components cannot be decoupled in the form that Eq.(2.9) represents the 
collision process. Decoupling can be justified, however, on the basis of the 
following argument. An examination of the terms on the right-hand side of 
Eq. (2.9) shows that the terms containing the M’s are production terms in 
which scattered molecules are produced at a rate equal to the product of the 
number density and the collision frequency or depending on 
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whether the collision process is self-collisions or ambient collisions. These 

/k 

terms are independent of f j and, thus, are decoupling. The terms containing 
f^ are attenuation terms in which the scattered molecules are attenuated as 
a result of subsequent collisions. Since the scattered molecules returning 
to the surface are traveling in the same direction as the ambient flow and in 
opposition to the outgassing flow, attenuation should occur primarily as a 
result of collisions with the outward flowing outgas molecules. Based on 
this argument, we drop the attenuation terms involving u^2 separate 

Eq.(2.9) into two parts containing contributions to the return flow due to self- 
collisions and ambient collisions: 

3fii afn 

+ V - fjj) (2.10a) 

d r 

„ Mj2-Pii£, 2 (2.10b) 

9t ^ r 

/\ A 

where is the contribution due to self-collisions and f is the contribution 
due to ambient collisions. For nearly free molecular flow, the attenuation 
terms should be small and, hence, any error due to their consideration should 
be small. 

Equations (2.10) are hyperbolic partial differential equations which can 
be integrated along characteristic lines in the direction of the velocity vector, 

V. These integrations will be performed separately in the following paragraphs 
for the contributions due to self-collisions and ambient collisions. 

Self-Collisions 

Integrating Eq. (2.10a) from a point at infinity, r , to a point on the out- 
gassing surface, yields the distribution function of return flow molecules 
on the outgassing surface; 
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* * 


f, , (v, ) = f 

11' w 


11 


(v, r^) exp 




( 2 . 11 ) 


At infinity the return flow aistribution function is zero, thus eliminating the 

first term in Eq. (2.11). The return flow rate q, at the surface is now 

°11 

obtained by integrating the return flow distribution function over all velo- 
cities directed back into the surface. 


q 


b 


11 



— * — k A — k — I 3 

V • e f , , (v, r ) d V 
n 11' w' 


( 2 . 12 ) 


where is an inward directed surface normal. Expressing the velocities 
in terms of absolute velocities, v, directions and the solid angle, uu , enclosing 
the directions yields: 


00 

■11 = //” 


j 1 ^ ^ 

V COS9 f, , (v, r ) dvdai 
11 ' w 


Zn 0 


(2.13) 


where (j) is the angle between the surface normal e^ and the velocity vector v. 
The solid angle integration is over the 2n steradians subtending the half- space 
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r 1 

j 

' 1 


outward from the surface. Combining Eqs. (2.11) and (2.13), along with the 
definitions of the Maxwellian in Eq.(2.3), yields: 



m 


2k0 


1 .2 

(V - Ujj) 

11 



V d V dr^ da.' 


(2.14) 


where the r' s are distances along the integral path. As noted earlier, atten- 
uation is small and, hence, the attenuation factor may be approximated by: 



Equation (2.14) now becomes; 


(2,15) 


00 / m \^/2 

11 °// (^^i) "• 


exp 


""I 2 2 

2k^ ^1 ‘"ll 


Zn 0 


00 

/ 


2 ’ ”"l 2 

- 2^j - ^n^osa^^) 


dv 


r °° 

/ ■'ll'*'"/ 


j 2 

''“P|- 2k5[7 ''' ■ "u '°“ll> 


dv 


dr'do> (2.16) 
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where ^ is the angle between the mean flow velocity vector Uj ^ and the 
velocity vector v. 

Integrating over the velocity yields: 


00 


dju 


(2.17) 


Ztt 0 


where is the directional distribution function (scattering pattern) of the 
scattered molecules, and gj ^ is an attenuation term. These parameters are 
defined by: 


11 


Ztt 


'^2 2 

, -u , sin a, , 

J- ^ ^ ^ * 

3/2“ ^ 


Uj J COSttj J G 


'^2 2 
-u j j cos j 


(1 + cos^aj j) 1 1 + erf (tij^ j cosajj)j 


(2.18) 


and 


/■'ll'’'" ~2 . 2 

0 “n 


■11 


27t^/^ ^2k0jj/mj 


~2 2 
Uj^ j cos a j j 


+ Ujj costtjj |l + erf (Ujj cosOj^ j)| 


where 


~ ”^1 
"ii ~ y 2k©jj "'ll 


(2.19) 


( 2 . 20 ) 
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The remaining integrations over r' and oj in Eq. (2.17) consist of 
simple volumetric integrals over the half-space outward from the outgassing 
surface (or the surface receiving the return flow). 

The parameters u^j^ and 0^^ depend on the local density and distribution 
function of the outgassing molecules (Eqs.(2.4)), and, so, will in general depend 
on the shape of the outgassing body and will vary with position in the flow field. 
The collision frequency I'jl assumed to be given by: 

•'ll = 

where is the outgassing surface temperature. In terms of the collision 
cross section and mean thermal velocity v^, this is equivalent to: 

1,111 VjOjjn^ (2,22) 

vthich coir*^ ares with 1.414 st Maxwellian gas (a gas in local 

tivermal equilibrium) and 0.806 v^ for a molecular beam. Our assumed 

collision frequency is thus seen to be a reasonable value, since the molecular 
flow would be expected to be much like a Maxwellian ^;as near the surface and 
a molecular beam far away from the surface. 

Ambient C Visions 


Equation (2.10b) is integrated in the same manner as in the preceding 
discussion to yield: 


12 


00 


Zn 0 


(2.23) 
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who re 


and 


. 2 

•“12 “ U 


12 - ^ 


Uj2Cosaj, 6 


~2 2 
■"12“°' “l2 


+ ^ (1 + 2uj^2COs^aj2) [l -I- erf (Uj 2 cosaj 2 ] 


(2.24) 


/ 


l/ljdr' 


1 0 


-2 .2^ 
-Ui2Sin aj2 


‘12 


27t^/^ ^2k0j2/i^l 


~2 2 

-Ui2C0Sj2 


>f7uj2Cosaj2 [: 


cosa,.~ 1 + erf (u,., cosa, 


(2.25) 


These equations are nearly identical to those derived previously for self- 
collisions. Note that the double subscripts ’* IT* are replaced by 12’’ every- 
where except in the collision rate pj^in Eq.(2.25). 

According to Eqs. (2.4), the parameters Uj 2 stnd 0j2 depend on both the 
outgassing flow properties, which are a function of position, and the ambient 
flow properties, which are essentially constant with position. As noted earlier, 
the relations reduce to the simple form given by Eq. (2.5) for the case of hyper- 
thermal flow. Since the ambient atmospheric molecules are traveling at very 
high velocities compared to the outgassing molecules for the orbiting space- 
craft we are considering, the use of these simpler relations in our analysis 
appears to be justified. 

The collision frequency 1/^2 again be assumed to take on the form 
assumed earlier in BGK model solutions: 
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U 


12 ' 


kT 


12 "2/^ 


(2.26) 


In this case, however, the temperature and viscosity must be characteristic 
of the collision process occurring between ambient and outgas molecules. 

A 

We will define the temperature T ^2 kinetic temperature characteristic 
of the energy of the flow (again neglecting the motion of the outgas molecules): 


A. 

T 


12 ■ 


”' l 2'^2 

3k 


(2.27) 


where ^ composite molecular mass for the two species. The viscosity 

11 is then defined in terms of this kinetic temperature and the collision cross 
section for collisions between the two species: 


M = 



(2.28) 


The collision frequency ihen becomes 


12 




12 ^ 2^2 


= 1,023 ^i2^2*^2 (2.29) 

Note that the composite mass ^^^2 bearing on the collision 

rate. The resulting collision frequency compares very closely with 1 ^ 22 "^ ^12^^2^2^ 
the relation for hyperthermal collisions between hard spheres. The BGK model 
approach is therefore seen to lead to physically reasonable results. 
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Section 3 

APPUCATION TO SPHERICAL SPACECRAFT 


The outgas contamination return flow model developed in Section 2 is 
applied in this section to a spherical spacecraft. The contributions due to 
self -collisions and ambient collisions are treated separately in the following 
paragraphs. 


3.1 SELF -COLLISIONS 


Consider the geometry depicted in Fig. 3-1. The free molecular dis- 
tribution function f ^ q outgas molecules at the point P is the Maxwellian 
distribution function corresponding to the sphere surface temperature T for 
velocity vectors enclosed by the solid angle Q subtended by the sphere. Out- 
side of this solid angle the distribution function is zero. The magnitude of the 
distribution function is determined by the outgas rate such that 





(3.1) 


on the sphere surface. The Maxwellian distribution func tion satisfying these 
requirements is 



= 0 


V not in Q 


(3.2) 


Using this distribution function, we obtain the following parameters; 


"1 


-fh,0 

n 
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(3.3) 


‘11 


57 /"‘i ,0 

1/2 


V 


n 


■ &) 


*^2 ~ a /-^2 T 

r - r ^ r - 1 


^^1 = “3F / 


= T 


w 


1 - 


. ~2 ~^\ 
l7r(r -r^r -1)J 


= T 


w 


- I ^11) 


(3.4) 


(3.5) 


I ^ 

where r = r/R, Uj^j^ = — '^ll ^ vector in the direction 

1 w 

outward along r. These parameters are all seen to be functions of the radial 
distance, r. When they are included in the scattering pattern function f^^ and 
the attenuation function Eq* (2-17) will yield the contribution to return flow 

due to self-scattering. Including these parameters in Eq. (2.17) and collecting 
into dimensionless groups yields: 

7t/2 


= Ztt Kn 


11 


■n! c. 00 

11 / sin0 cos^ / ^11 


dr' d(^ 
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7T/ ^ 00 

- Ztt Kn J* sin</)cos0 ^ ®11 




‘11 “ m, R 




= l-rr.l/' 

®11 =®llAw 


f ^ 


■?i ^ 


u, , sin a, , I ~2 
11 11 


+ ^7 UjjCosttjj |l + erf (U|^ j cosaj^ j)j 


The parameter Krij j may be considered a Knudsen number character- 
istic of self-collisions. The parameter n^ may be considered a dimension- 
less number density, and the integral 




may be considered a dimensionless column density. This integral may be 

• • 

carried out in closed form by approximating n^ by 
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nj = 1 




(3.9) 


This approximation was made with good results in the study document d in 
Ref. 6. The resulting expression for N is 


N 



/ r ' + cos^\ 

y sin<^| j 

+ 0 - 

7t/2 



2 , 

- 


r ' + cos<^ 

sin 0 

+ 1 

- cos(j> (sin^i^ + ^ 



r 

r 

J 



I 


tan 


-1 


( r ' + cos^ I 
sind / 


+ ^ - 7t/2 



The radial distance r can be expressed in terms of r ' by 


1+2 cos</> r ' + r 


(3.10) 


(3.11) 


The remaining integrals to be performed in Eq.(3.6) are simply a volu- 
metric integral in axisyxnmetric spheric coordinates. These integrals, v/hen 
evaluated numerically, yield the final expression for return flux due to self- 
collisions: 


= (0.01554 - 0.004058 Kn'})Knj} 


(3.12) 


This expression is plotted in Fig. 3-2 as a function of the inverse Knudsen 
number, Kn^ and coii.percL with results obtained from a simple test particle 
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Monte Carlo program. This Monte Carlo program was developed solely to 
provide a source of comparison for the BGK model results. The BGK model 
results are seen to be somewhat less than the Monte Carlo results, buc still 
well within the same order of magnitude. 


3.2 AMBIENT COLLISIONS 


The geometry involved in the integration of the BGK equation for ambient 
collisions is represented in Fig. 3-3. Spherical symmetry is no longer applic- 
able in this case because of the directed flow of the ambient freestream. The 
return flux at a point Q on the sphere surface now depends the angular dis- 
placement </>' away from the forward position (stagnation point), “^he collision 
rate Eq.(2.29), and the local number density n^, given by Eq.(3.3), 

may be included in Eq. (2.23) to yield the dimensionless form: 


= 1 


12 




^’^ 22^22 


00 


y*cosc^.fi2 J 


d r' du: 


27T 


00 


-■yjir Kn”J J coF<t> f 


^ du; 


(3.13) 


2 ?: 0 


where 


^2 V nij 


(3.14) 


Kn 


22 


V ^>^2 ”22 ^ 


(3.15) 


and 
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Fig. 3-3 - Geometrical Representation of Scattering by Ambient 
Atmospheric Molecules 
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~2 2 
-u^2COSj2 


+ V^ '^12^°®*^12 ^ erf (Uj^2 


(3.16) 


The parameter is the speed ratio, and Kn 22 is the freestream Knuden 
number. Using the values Uj 2 a^nd 0^2 given by Eqs.(2.6), the functions fj 2 
eind gj 2 become: 


'12 


2n 


3 . 2 

1 ■ I “l2 


3 2 

3 2^°® “l2 

'2' cosa^2 G 


VI 




(1+3 cos ttj 2 ) I ^ IV 7 I 


(VI 


12 . 


rr (1 + "^1^2) ~ , 

®12"12 S, 


12 


(3.17) 


(3.18) 


where 


3 . 2 

- 2 sin a ^2 


■12- ^ 


3 2 

I 


cos a,^ 


+ cosa ^2 [l + erf (^ cosaj 2 )] 


(3.19) 


After inclusion of the parameters defined by Eqs.(3.18) and (3.19), Eq. (3.13) 
becomes: 
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(3,20) 


The integration over the radial distance r * in the first term of Eq. (3.13) 
is just the uimensionless column density at infinity which from Eq.(3.10) is: 


00 


/ 


nj d r' 


M - 1 « ^ 4 . 1. ^ 1 cos4> 2 

~ 2 I sin0 8 . 5 . " 4 . 2 1 8 

\ s’.n 9 sin 9 


cos9 \ 
sin^9 / 


(3.21) 


The integral over r' in the second term is: 


00 


/ 


n, N dr' = 4 N 
1 2 00 


(3.22) 


Equation (3.20) now becomes; 


q, = 1.023 Knll 

^12 °22 


/ 


S, I COS0 f 1 o N duu 
by 12 00 

Ztt 


- I + m/rr2)Kn;} cos9h^2N^ cL; 

Ztt 


(3.23) 
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The angle ct-^^ given by (see Fig. 3-3): 


-1 


a ^2 = cos” (cos(/>' cos^ - sin0' sin^cosG) 


3.24) 


At the forward station = 0), 0:^2 = ^ and the integral over the Ztt steradians 
in Eq.(3.23) becomes a single integral over ^ from 0 to tt/Z. At other stations, 
the integral requires a double integration of ^ from 0 to tt/Z and 0 from 0 to 
Ztt. In all cases, the integration must be performed numerically. The result- 
ing form of the solution is: 



^ Kn”2 Sj^F(<^') - (1 + m^/m2) Kn”}G(«')| 


(3.25) 


where the factors F and G depend on sphere station. The computed values of 
F and G are given in Table 3-1 compared to results obtained previously using 
a ’’first collision” model approach (Ref. 6). The ’’first collision” res ults did 
not consider an attenuation factor, G* The BGK and ’’first collision” model 
results are seen to compare fairly closely over the entire range of sphere 
stations with the exception of the rear stations. The first collision model 
does not allow upstream scattering, and, hence, the return flux approaches 
zero at the rear station (180 deg). The BGK model, on the other hand, impli- 
citly includes a combination of first and multiple collisions which allow some 
upstream scattering. The BGK model thi^s results in a nonzero return flux 
at the rea station. 
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Table 3-1 

COMPUTED PARAMETERS FOR BACKSCATTER TO SPHERE 
DUE TO COLLISIONS WITH AMBIENT ATMOSPHERE 


(d-g) 

F 

(BCK) 

G 

(BGK) 

F 

("First Collisions,") 
Ref. 6 

0 

,607 

.304 

.56'^ 1 

10 

* 606 

.304 


20 

.587 

.298 

.534 

30 

,554 

.287 

.501 

40 

.508 

.268 

.459 

50 

.451 

.244 

.410 

60 

.385 

.216 

.358 

70 

.315 

,184 

.303 

80 

.245 

.150 

.248 

90 

.181 

.118 

.196 

100 

.128 

.0890 

.150 

no 1 

.0877 i 

.0659 

.110 

120 

.0590 

.0485 

.0757 

130 

.0399 

.0361 

.0479 

140 

.0278 

.0278 

.0267 

150 

.0205 

.0223 

.0122 

160 

.0162 

.0190 

.0039 

170 

.0140 

.0172 

.0005 

180 



.0134 

.0167 

.0000 
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Section 4 

APPUCATION TO FLAT PIATE DISK 


The return flow model developed in Section 2 was also applied to a flat 
plate disk oriented normal to the direction of the ambient flow. In this appli- 
cation the return flow was calculated at points at various elevations above the 
surface of the disk. The method of application was similar to that described 
for the sphere, but some liberties were taken to simplify the analysis. Con- 
sider the disk geometry depicted in Fig. 4-1. The return flux is to be deter- 
mined at the point Q located at an elevation, H, above the center of the disk 
surface. The free molecular distribution function, f^ q, for the outgas mole- 
cules at the point P is the same as given for the sphere in Eq. (3.2) for velo- 
city vectors enclosed by the solid angle subtended by the disk- Outside of this 
solid angle, f ^ q is zero. The outgas flow parameters are easily obtained by 
integration of the distribution function along the line perpendicular to the disk 
center (</>* = 0). For points off this line, however, the integrations become much 
more complicated and are probably not obtainable in closed form. To simplify 
matters it was assumed that, for points within a distance of one disk radius of 
the disk center (r < R), the local number density, n j, is equal to the centerline 
density at the same distance from the disk center point. For points beyond 
one disk radius (r >R), the density was calculcated based on a far field approx- 
imation (r »R). In the far field, the outgassing disk is approximated by a point 
area source emitting with a Lambertian (cosine) distribution. The mean flow 
velocity u^^ and temperature 0^^ are not strongly dependent on position, and 
so are assumed to always take on centerline values. The resulting outgas flow 
parameters are: 





1/2 


CQS<(>' 



r > R (4.1) 
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Fig. 4-1 - Geometry Involved in Back-Scattering to Points On and Above 
Flat Plate Disk 
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1/2 


~2 

r - r 


VF7 


(4.2) 


1 + 1 


0,, = T 
11 w 


' -If 


~2 

(r - r ^ r 


+ 1+1 


T (1-1??, 
w ' 3 11 


(4-3) 


where r 


r/R 


•“u =■^1 


m, 


u. , and e is a unit vector in the direction 
1 i r 


2kT 

w 

outward along r. The integration of the BGK model equation for the flat plate 
disk follows closely the process described in Section 3 for the sphere. The 
resulting expressions ijr the self-scattering and ambient scattering contri- 
butions to the outgas return flux are: 


% 


11 


F(H) - G(H)KnJ^ 


Kn 


-1 

11 


(4.4) 


12 


''12 ^ -1 
— — Kn_^ 
0^2 22 


-1 


S.^F(H) - (1 +m^/m2)Kn j j G(H) 


(4.5) 


The parameters F and G are plotted for the case of self-collisions in Fig. 4-2 
and for the case of ambient collisions in Fig. 4-3. 
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Fig.4- 


2 - Return Flux to Points On and Above Flat Plate 
Disk Due to Self-Collisions 
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Section 5 

APPLICATION TO SPACE SHUTTLE 

The equations developed in Section 3 for a sphere and Section 4 for a 

flat plate disk were used to estimate the backscatter of outgassing from the 

Space Shuttle vehicle. The Space Shuttle Orbiter is approximately 35m long 

and 24m wide at the base. We assumed an overall equivalent diaineter of 

20m for use in the sphere and disk return flow equations. Calculations were 

made based on a circular orbit altitude of 400 km. A surface temperature of 
o 8 2 

300 K, an outgas rate q^ of 10* g/cm /sec, and an outgas product molecular 

weight of 100 were assumed. Atmospheric data were taken from Ref. 7 which 

_ 26 2 

provided a value of 41.9 10* cm for the ambient c\oss section Based 

- 16 2 

on molecular data in Ref. 8, a value of 43.0 x 10 cm was assumed for the 
collision cross section, of the outgas molecules at 300^K, This value was 

also taken as the cross section collisions between outgas and ambient 

molecules at energies corresponding to a reference temperature T of 300 K. 
The cross section corresponding to orbital interaction velocities was esti* 
mated based on the temperature variation of viscosity: 


■(" • ?) 


(5.1) 


/A 

where T ^2 kinetic temperature defined by Eq. (2.27) and n is the ex- 

ponent in the power law variation of viscosity with temperature. Assuming 
a composite molecular w'eight of 50 for the interacting species, and a value 

/S. 

of 0.8 for the exponent n (Ref. 8), the interaction temperature T ^2 found 
to be about 100,000^K, and the cross section ratio is about 0.2. 


The self-scattering return flux was calculated and found to be 2.25 x 
10* g/cm^/ sec (1.35 x 10^^ molecule s/cm^/ sec) using the sphere formula 
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- 12 2 10 2 
and 3.80 x lO” g/cm /sec (2.29 x 10 molecules/cm /see) using the disc 

formula for the zero elevation. The self-scattering contribution is the same 

for all stations on the sphere. The ambient scattering return fluxes on the 

sphere were found to be 4.83 x lO”^^ g/cm^/sec (2.90 x 10^^ molecules/cm^/ 

- 12 2 9 

sec) at the forward station (<i>' = 0) and 1.07 x 10” g/cm /sec (6.40 x lO^mole- 
2 

cules/cm /sec) at the rear station (<^»' = 180 deg). Note that the predicted 
ambient scattering return flux is about 20 times greater than the self- scattering 
return flux at the forward station and about half the self- scattering ret.irn flux 
at the rear station. The fact that the predicted ambient scattering return flux 
at the rear station is within the same order of magnitude as the self- scattering 
return flux is noteworthy, since previous studies based on first collision theory 
did not allow any upstream scattering in ambient collisions. The validity of 
these upstream scatt«. .ing predictions wi’. be discussed later in Section 6. 

The ambient-scattering return flux on the disk was found to be 7.43 x 
10 g/cm^/sec (4.46 x 10^^ molecules/cm^/sec). The predicted return flux 
to the disk is seen to be about 50% greater than that to the sphere at the forward 
station . 

The contamination control criteria given in Ref. 9 specifies a return flux 
less than 10 molecules/cm /sec. The predicted return fluxes for the 400km 
orbit satisfy these criteria. 
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Section 6 

PHYSICAL INTERPRETATION OF THE BGK MODEL SOLUTION 

In light of the BGK model development described in Section 2, the phys- 
ical interpretation of the BGK model solutior. can be stated as follows. Scattered 
molecults are produced locally at '^ll points in the flow field at a rate equal to 
the local collision frequency. These molecules take on velocities and directions 
consistent with a gas in equilibrium at a hypothetical temperature and with a 
superimposed mean flow velocity. The hypothetical temperature and mean 
flow velocity are determined such that overall momentum and energy are con- 
served in the collision process. The BGK model scattering pattern, therefore, 
satisfies the basic requirements of mechanics, but it is not unique in doing so. 

The scattering pattern for self-collisions for cutgassing from a sphere 
is shown in Fig. 6-1, with distance away from the sphere center as a parameter. 
The scattering is seen to be highly directed in the direction of the mean flow 
away from the sphere, with the distributions becoming more directed as the 
distance increases. The scattering pattern approaches a constant distribution 
in the far field. Although the distributions are highly directed in the direction 
of the mean flow, there is in all cases a small, but nonzero, component of back- 
scattering. Mathematically, this is brought about simply because the high 
velocity "tail" of the assumed Maxwellian distribution function contains some 
molecules with upstream components greater than the superimposed mean flow 
vciucity. There is a good deal of physical basis ft r the predicted upstream 
scattering for the case of self-scattering from a sphere, since even "first" 
collisions of molecules emitted at high angles from the surface normal can 
cause deflections back into the surface. 

More difficult to explain is the predicted upstream scattering for the 
case of ambient collisions. The predicted BGK scattering pattern for this 
case is shown in Fig. 6-2 compf red to the cosine distribution characteristic 
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Fractional Scattering Distribution per Unit Solid Angle, f 


Fig. 6-1 - BGK Model Scattering Pattern for Return Flux 
Due to Self -Scattering 
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Fractional Scattering Distribution per Unit Solid Angle, f 

Fig. 6-2 - Scattering Patterns for Return Flux Due to Collisions 

with Ambient Atmosphere — BGK Model vs Cosine Law 
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of ”first^’ collisions in hyperthermal flow. Again, there is a small, but non- 
zero, component in the BGK predicted upstream scattering. Since first*’ 
collisions between outgas molecules and freestream molecules do not permit 
upstream scattering, the physical interpretation of the upstream scattering 
components is that the overall collision process consists of a spectrum of 
collision events including not only ’’first'* collisions, but also more compli- 
cated multiple collisions. Based on a cosine distribution, a small fraction 
of molecules will be directed in a near-lateral direction with respect to the 
ambient velocity vector after a ’’first” encounter. Some of these molecules 
will be directed in opposing directions such that ’’second” collisions can occur 
between these molecules to produce upstream as well as downstream deflec- 
tions. Other more complicated collision events can also produce upstream 
deflections. Since the BGK model made no attempt to incorporate a rational 
combination of these collision events in the collision process, and since the 
upstream components form a relatively small portion of the overall scattering 
pattern, it is probably safe to say that there is a considerable margin of error 
in the BGK model prediction of upstream scattering for the case of ambient 
collisions. Probably the truth lies somewhere between the BGK and cosine 
distributions, which means that the actual upstream scattering for ambient 
collisions is less than the BGK model prediction. 
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Section 7 

CONCLUDING REMARKS 

The BGK model has been shown to be a useful device for estimating 
return flux for both self and ambient- scattering. The functional form cf 
the return flow equations developed from the BGK model has the same form 
as equations developed previously by Scialdone (Ref. 10) and Robertson (Ref. 6) 
using the first collision model appt.oach. The present approach, however, in- 
cludes an attenuation term and provides for a small amount of upstream scatter- 
ing in ambient collisions. As discussed in Section 6, the upstream scattering 
in ambient collisions is the primary area of doubt in the BGK approach. It was 
shown in Section 6 that a mechanism for upstream scattering exists, but at this 
point, it is impossible to quantify it. It is certainly likely that the BGK esti- 
mates provide an upper bounds. The return flux calculations in Section 3. 
therefore, may be considered reasonably accurate for self-collisions and fcr 
ambient collisions on the forward side of the sphere. The wake side calcula- 
tions for ambient scattering, however, must be regarded as upper bounds only. 

A significant point was noted by Naumann (Ref. 11) that the return flux 
of an outgassing species is dependent on the total molecular outflow. The 
flux q in the Knudsen number Kn. , for self-collisions , therefore, would 
be the total outflow rate. Thus, a contaminant species diluted in a large 
efflux of benign species would be returned at a higher rate. 
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